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The leading (Gaussian) fluctuation correction to the weak coupling zero temperature BCS super- 
conducting gap equation is computed. We find that the dominant contribution comes from the high 
energies and momenta (compared to the gap) and gives a correction smaller by the weak-coupling 
factor gNo than the mean-field terms. This correction is small due to cancellation of singular 
contributions from the amplitude and phase mode at high energies and momenta. 
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I. INTRODUCTION 

The Bardeen-Cooper-Schreiffcr (BCS) theory of 
superconductivity^ is important both as an empirically 
highly accurate and successful model of a nontrivial phys- 
ical phenomenon and as a paradigm for theoretical study 
of a wide class of models involving a logarithmically di- 
vergent susceptibility. This paper focusses on the latter 
aspect. From this point of view the essence of BCS the- 
ory is the observation that in a generic, weakly coupled 
fermion system in d > 2 dimensions, all susceptibilities 
arc non-negative and remain finite as temperature T — > 
except for particle-particle susceptibilities such as 



Tr 



(1) 



^^(T,r)V'T(T,r),^+(0,0)7A+(0,0) 



which diverge logarithmically as q, v and temperature T 
tend to 0. (By the phrase 'generic' we mean to rule out 
for example the nesting and van Hove instabilities occur- 
ring in particular models at particular band fillings). 

BCS showed that the logarithmic divergence of x sig- 
nalled the appearance, at a transition temperature Tc, 
of a new order parameter A ~< V'J.V'T >• They further 
argued that Tc, the magnitude of A and many other phys- 
ical consequences of the ordering could be accurately de- 
termined by mean field theory. The stunning agreement 
between predictions of the BCS theory and data on con- 
ventional superconductors lends very strong support to 
this view. 

The approach pioneered by BCS has been applied by 
many workers in many contexts involving logarithmically 
or more strongly diverging susceptibilities; for example to 
spin and charge density wave instabilities driven by nest- 
ing effects in half filled particle-hole symmetric bands^, 
or, in the two dimensional case, at fillings that give rise 
to a van-Hove singularity in the density of states'^. 

Some aspects of fluctuations have been understood in 
detail. Around the time BCS theory was developed it was 



understood that in d < 4 dimensions and for tempera- 
tures sufficiently near to Tc, nonlinear interactions among 
very long wavelength order parameter fluctuations would 
invalidate a mean field treatment of the thermally driven 
transition, and the (rather modest) temperature window 
within which conventional superconductors would exhibit 
non-mean-field behavior was estimated^-. The effects of 
order parameter fluctuations with momenta and energies 
smaller than the inverse correlation length and A, re- 
spectively, on the near- Tc properties of conventional su- 
perconductors in different dimensions has been examined 
in detail^ and found, for example, to give a small negative 
correction to the mean-field value of the order parameter 
A and in d = 2 to have a larger effect on Tc- Engelbrecht 
and co-workers used a functional integral method to in- 
vestigate the effect of long wavelength fiuctuations on the 
BCS-Bose-Einstein crossover—. Varlamov et al— studied 
the effect of these long wavelength fluctuations on the 
normal state of layered superconductors in the context 
of high-temperature superconductivity. 

In this paper we investigate corrections to the BCS 
mean field approximation due to Gaussian fluctuations. 
We use a standard functional integral formalism to study 
the fiuctuations over a wide range of momenta and ener- 
gies at T = 0, where the order parameter is believed to 
be well developed, and thus the fluctuations corrections 
are supposed to be small. Suprisingly, this basic question 
seems not to have been addressed in the literature. We 
calculate the corrections and flnd that they are, indeed, 
small compared to the mean-field terms. There are three 
important features: 1. The dominant correction comes 
from energies and momenta high compared to the mean- 
field gap. 2. The dominant contribution to the correction 
comes from processes in which the electrons are scattered 
nearly parallel to the Fermi surface. 3. The fluctuation 
corrections from each of the amplitude and the phase 
mode diverge logarithmically at high energies, but with 
opposite sign so these divergences cancel leaving a small 
overall correction. 

The rest of this paper is organized as follows: in Sec- 
tion II we present the formalism to be used. In Section 
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III, we evaluate the polarization kernel and its deriva- 
tives, and use it in Section IV to calculate the fluctuation 
correction. Section V is a discussion and conclusion. 



II. GENERAL FORMULAE 

We consider for definiteness a model of fcrmions in 
d > 2 spatial dimensions with energy dispersion Sp — 

/2m — /i (as will become evident below, our results 
may be trivially generalized to more realistic dispersions 
provided nesting and van Hove singularities are absent). 
We take the fermions to interact via a short ranged in- 
stantaneous attractive interaction parametrized by a co- 
efficient 5 > 0. We follow ShankarS and consider only 
states within a cutoff A << pF of the fermi surface, ex- 
pand the density of states A^(e) = / L (5(£ — £p) around 



the Fermi energy as 



N{€) = No + Nie/vFl^ + 



(2) 



where Nq and A'^i are of the same order of magnitude. 
The Hamiltonian thus becomes 



H = J depEpC^^Cpa 



(3) 



-g I {dpi...dp4)c+^^.c+^^Cp.,iCp^'i 



where (dp) — ^J!^d f^nd the prime on the integral indi- 
cates that all momenta are restricted to the shell —A < 
\p\ — Pf < ^ and that there is a momentum conserving 
delta function. 

To analyse the model we write it as a functional 
integral^, decouple the interaction via a Hubbard- 
Stratonovich transformation and perform the integral 
over the fermionic fields obtaining for the partition func- 
tion 



J VA*{T,r)VA{T,r)e' 



with action S given by 



S ^ Trln(-9^ - H{A{t, r))) - J drdr 

0V 



|A(r,r)p 



(4) 



(5) 



and 



i/(A(r,r)) ^ depEpC+^Cp, (6) 

(A(r, r)2P+ (r, t) V'j*' (r, t) + H.c 



drdr 



13V 



The BCS mean field theory corresponds to a saddle-point 
approximation to Eq^ at the saddle point we have 



A(r,r) = A* (T,r) = Ao, 



(7) 



and the saddle-point approximation to the action is 

5o(Ao)-5o(Ao = 0) = /5y ( NoAlln^^ - ^ + . 

(8) 

where the ellipsis denotes terms of the order of Aq with 
coefficient of order unity, and vp = dsp/dp at pp is the 
Fermi velocity. Our calculations are valid in the weak- 
coupling regime gNo << 1, and these omitted terms are 
small by at least one power of gNg relative to terms which 
we retain. Finally, the saddle point value of Aq is fixed 
by extremizing S with respect to Ag yielding the familiar 
BCS gap equation 



9 \ Aa 



(9) 



where again the ellipsis indicates terms of order unity. 
Thus, as is well known, within the BCS approximation 
the T = gap value Aq is determined only to logarithmic 
accuracy, i.e. ln(A/Ao) is known up to terms of relative 
order unity. 

To study fluctuation corrections we write: 



A(T,r) = Ao + ?7(r,r). 



(10) 



We write the fluctuation rj in terms of its real and imag- 
inary parts 



v{T,r) = 771 (T,r) +iT]2{T,r), 



(11) 



which, for the gauge in which Aq is real, correspond to 
the amplitude and phase fluctuation respectively. 

Substitution of Ea llUl into Eq^and expansion in pow- 
ers of 5 yields Z = J Drj*{T, r)Dri{T, r)e'^o+^S ^-^^^i 



vli-i'^, q)vb{i'^, q) + 

where the ellipsis denotes higher powers of rj, a,b — 1.2 
and 



{-iY+^ r duj r /dk 
2 J 2^ J [2^^ 

Tr[G{iui+,k+)TaG{iuj-,k^)n] 



(13) 



Here lo± = lu ± ^ and k± = k ± ^ while G is the su- 
perconducting Green function in Nambu matrix notation 
and Ta,a = 1,2 are the Pauli matrices in particle- hole 
space. Note that in the instantaneous interaction model 
it suffices to place a cutoff on the momentum integral; no 
frequency cutoff is required. 

We see from Eq^jthat the propagator corresponding 
to 77 is of order g (except at long wavelengths and low 
energy) so an evaluation of Z via an expansion in powers 
of 1] yields a series expansion in powers of gN^ for the 
free energy and other physical quantities. For example. 
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restricting the expansion to Gaussian order leads to a 
correction to the free energy given by 



Fi = -TS, 



(14) 



^^E/ (§) Trln(l + gn(z^.,g)) 

A ^ 

^^E/(|) ^ln((l + ^;no)^ (15) 



Here 



n EEE Hoi + HlTl + n2T2 + II3T3 



(16) 



is the polarization matrix written in terms of the Pauli 
matrices. The correction to the saddle point equation 
arising from the Gaussian fluctuation term is 



and Hi = 0. 

First, we analyze 113. We linearize the fermion disper- 
sion around the Fermi surface 



(22) 



where is the Fermi velocity, and ^ is the cosine of the 
angle between k and q. For a spherical Fermi surface in 
d dimensions, vp is a constant, and 



N{e)de 

-vfA 



with 



Then 



Kd = 



2r(f) 



(23) 



(24) 



(17) 



(1 + .9no)g^ - g^(ni|g» + + n3^) 
(i + .gno)2-.g2(n2 + n2 + ni) 

In the next section, we shall estimate the relative contri- 
butions of Eg 1171 and EqlHlto the gap equation. 



III. CALCULATION OF n 

This section computes the polarizibilities appearing in 
Eq 1171 To simplify notation we work with the action 
density (i.e. divide by (3V) and omit the subscript of 
Ao. ^ 

Using the explicit mean- field form for G, namely 



G{iuj,k) 
we find (here e± = ek±q/2) 



tj2 + e'^(^k) + A2 



(18) 



Tin 



did f f dk\^ 



2tt J \2tt J 

cj+aj_ -I- e+e_ 



(19) 



Ho 



duj f f dk\^ 
2^ J [2^) 



(20) 



4 + 

A 



A2] 



„ , duj f I dk 

= '2^ (2^ 



(21) 



A2 



H + 4 



A2] 



Ilsiiiy, q) = NoJ ML_Z!l^/(^, A; A) (25) 



with 

(26) 

-vpA 

A2 



I{ry,VFqtJ-,A;A) ^ J ^ J de 



+ e+2 _^ A2] [w_2 + e_2 + A2] 

If << vpA and q << A then I depends on i^, q only via 
the combination 



2A 



(27) 



so we may evaluate the e and lo integrals at q = and 
v — 2Ar obtaining 



/(r;A) EE 



vp-A 00 

de 

J 2tt 

-vpA —00 



(28) 



{uj + Arf + + {uj ~ Ar 



1 



2rVr 



Infr + VT 



(29) 



where the last approximation applies for rA << vpA . 

Next, to analyze Hq, we separate no(0, 0) from the 
formula by adding and subtracting ^ [1^+'^ -I- -f- A2] -|- 
i [uj'i + £2^ + A2] to get 



Uoiiiy,q) = Ho (0,0) 



(30) 



-No 1 '^^'-^'^' (2.^ + l)/(.,A) 
Hd 
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For the static uniform polarization, we obtain for A << 



Ho (0,0) = -iVo 



2^ 



de- 



(31) 



Nn In 



(32) 



where again the ellipsis denotes terms of order unity. 

Finally, we analyze 112. We see that U2{ii',q) 
as 1/ ^ 0, so 



n. 



—ivJiiu^ q). 



(33) 



To the accuracy with which we discussed II3 and Hq, we 
may say that J is a function of r, and we evaluate it at 
q = 0,v = 2Ar. Thus, 



J 



vpA 

de 
N{e)e 



dio 
2^ 



(34) 



[{uj + Ar)2 + £2 + A2] [{uj - Ar)2 + + A^] ' 

J vanishes in case of particle-hole symmetry, so the con- 
stant term in the density of states (Eq ^ does not con- 
tribute. The Lu integral can be done analytically and we 
find 



J 



d/x(l — /i2) 2" 
Kd 



(35) 



vpA 



de- 



-vpA 



Ve2 + A2(e2 + A2 A2r2) 



We separate out the logarithmically divergent term 
and perform the integral over e obtaining 



J = 



2vfA 



rfM(l-M')' 

Kd 



(36) 



In 



The derivatives dH/dA^ needed for Eq El n^^^y be 
computed by straightforward differentiation (recall that 
r = \J ^ (vpq[L)^ jlA). We summarize the results: 



no = -A^olnf^) 



V A 



(37) 



d/i(l - p?)-^ (2r2 + 1) ln(r + ^.fF^^) 



Hd 



Ho 



r - In -7— 



(38) 



dfi{l - fi'^)^ %/r2 + 1 ln(r + Vr^ + 1) 



Hd 



dUp 
dA? 



dfi{l - /i2)T^ ln(r + \/r2 + 1) 



2rVr2 + 1 



No rd^l{l~^?Y 



4A2 



Hd 



(39) 
(40) 



ln(\/r2+T + r) 1 



r(r2 -I- 1)3/2 r2-|-l 



dA2 
dA2 



4A2 



Hd 



4A2 J Krf 





(41) 
(42) 



(2r2 + 1) ln(V?'2 + 1 + r) 1 



r(r2 1)3/2 



+ 1 



IV. EVALUATION OF FLUCTUATION 
CORRECTION. 

This section uses the results of the previous section to 
evaluate EqEI There are two regimes in Eq (|17|l : small 
r {i/'^, {vfqY ^ ^-nd large r. We consider them in 
turn. 

(i) Small r: in this limit both Il2dIl2/dA^ in the nu- 
merator and in the denominator of H17|l are suppressed 
by the small factor (ly/vpA)'^ relative to the other terms, 
so we shall neglect the contribution from 112. In d = 2 
with the parabolic dispersion, iVi = identically, so 
this term vanishes altogether. Then the leading behavior 
of both 1 -I- gllo and gll^ is gNo/2; the leading behav- 
ior of c?no/dA2 is No/2A'^, and the leading behavior of 
dUs/dA'^ is reduced compared to Nq/2A^ by r"^ . Thus, 
the numerator behaves as (gA^o)^/4A2, and we factor the 
denominator as 

(1 + .gHo + 5n3)(l + glio - .gHg) ~ [gNpfl x (43) 

identifying the two factors with amplitude and phase. 
That way, we find 



(44) 



dSi 


A'^-^ fdv j 


Sd- 




1 


dA'^ \small 


Vf J TT J 






^d 






Sd-iA''^' 
d-K'^+'^vj, 


-No 


2'^ / 
^ 1^ 


^) 

vpAj 



with u = ^ifS- and v — since, for a spherical Fermi 



2A 



5 



surface, 



Sd-iA 



d-l 



(45) 



(ii) Large r: The leading behavior of the individual 
terms in the numerator and denominator of (|17|l is 



(l + .gno)5^ 



■lnr(46) 



In the first region, we neglect the u dependence of r, so 
the fj, (angular) integration becomes trivial, and we get 



d-i 



dv 



d-i 



K 2A 



1 



du 



No 



(54) 



1 7r(d-l)lni^' 



In the second region, we find, with logarithmic precision, 



Kd r2 



1 

/ 



Kd 



■ In 2-\/ + (u/i)^ = Inu, 



(55) 



and 



, dU2 



2vfA 



In 



dfiii-^i^y 



Inr 



Kd 



dng ^ {gNoY / f dfi{l-fi^)^ Inr 
dA2 4A2 I ./ r-2 



(47) 



(48) 



so 



1 



A'^-^Sd-i 



2A 



1 



+ (u/i)-^ KdUW 



1 



KdUU Inu 



2vfn'^+^ {d - l)Kd V 2A 



(56) 



(57) 



7r(d - l)Kd 



{l + gnoY - (5iVo)^ f / ML^f!)ZinH (49) 



(5n2)2 - -{gN,y 



(gn^f = (giVo)^ 



2vfA 



Kd 
2 

In' 



2 VpA. 



dn{l — lY) 2 In) 



(50) 



(51) 



We see that the leading contribution to both the nu- 
merator and denominator comes from Ho, so 



dSi 
dA^ 



Ai 



dv I I dq\ g-^ 



2-n J \2'K J 1 + gHo 
A'^-'^Sd-i f'd-i. 



du / dv 



(52) 
(53) 



1 , 9 

r d^ijl-fj.^)^ 1 



1 d-3 





where the prime means that we are integrating over large 
momenta and frequencies compared to the gap. Now we 
have to consider separately the regions u < v and u > v. 



so it is the dominant fluctuation contribution to the fluc- 
tuation correction of the gap equation. We see that it is 
by factor f;7Vo or l/ln('yFA/A) smaller than the mean- 
field terms, and that it is negative, so it decreases the 
value of the gap. This contribution will change only the 
prefactor in the solution of the gap equation. Other ef- 
fects, not considered in the present calculation, will also 
make corrections of the same order—. Comparing for- 
mula H54I) to formula H57II . we see that the contribution 
from processes that scatter electrons along the Fermi sur- 
face dominate other contributions by factor ln(u^A/A). 

It is interesting to note that the result (|57|l is small 
because of cancellation of large terms. Indeed, returning 
to p4|) . we can write the fluctuation correction to the 
gap equation ((T7jl as 



dSi_ 
dAf 



-(3V 



dv 
2^ 



Tr 



dU 

'dAf 



(i + ^n)- 



(58) 

At high frequencies and momenta, the leading behavior 
of the two eigenvalues of the matrix 



dU^ , 



(59) 



that is, of the amplitude and phase mode, is ±/ where 



1 + gn 
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A 



(60) 



SUMMARY AND CONCLUSIONS 



1 „ d-3 

r dfi{l-fi ) --1 Inr 



1 d-3 



llli 



Calculations similar to those leading to eqs H54|l - H57(l 
show that, with logarithmic accuracy, 



1 A 



d-l 



vfA/2A 
Sd-1 f _d-l 



d+1 



In V 

du I dv-^ 
mu 

d-l 



1 A-^-i Sd-i 1 1 (vF^y-' \n' ^ 



vi TT-^+i Kd 2 d - 1 V 2A 



In 



vpK 



Nn 



7r(d - l)Kd 



In 



(61) 



Thus each of the amplitude and phase mode gives a 
correction apparently large enough to call the BCS ap- 
proximation into question, but the two terms cancel in 
the trace. 



We have used a functional integral formulation to 
study fluctuation corrections to the weak coupling BCS 
mean field expression for the superconducting gap am- 
plitude. We find that in both two and three dimensions 
and in the weak coupling limit, the correction is smaller 
by the weak-coupling factor qNq or l/ln(w^A/A) than 
the mean-field terms. This correction comes from ener- 
gies and momenta large compared to the mean-field gap, 
and, more specifically, from the region in the phase space 
where the electrons are scattered nearly parallel to the 
Fermi surface. We note that individually, both the con- 
tribution from the amplitude mode and the contribution 
from the phase mode diverge logarithmically in this re- 
gion, but the divergencies have opposite sign, and thus 
cancel when we take the trace. 



Acknowledgements: A. J. M was supported by NSF- 
DMR-0431350. S. K. acknowledges the LANL DR 
Project 200153 and the Department of Energy, under 
contract W-7405-ENG-36, and A.I.L. NSF 0439026. We 
thank L. B. loffe for helpful discussions. 



^ J. Bardeen, L.N. Cooper, and J.R. Schrieffer, Phys. Rev. 

113, 982 (1959) 
^ A.W. Overhauser, Phys. Rev. 128, 1437 (1962) 
^ T.M. Rice and G.K. Scott, Phys. Rev. Lett. 35, 120 (1975) 
* A.P Levanyuk, Sov. Phys. JETP 36 571 (1959) 

V.L. Ginsburg, Sov. Phys.-Solid State 2, 1824 (I960) 
® L.G. Aslamasov and A.I. Larkin, Phys. Letters 26 A, 238 

(1968) 

K. Maki, Prog. Theor. Phys. (Kyoto) 39, 897 (1968); 40, 
193 (1968) 



R. S. Thompson, Phys. Rev. B 1, 327 (1970) 
® A. A. Varlamov, G. Balestrino, E. Milani, D.V. Livanov, 
Advances in Physics, 48, 6, 1 (1999) 

J.R. Engelbrecht, M. Randeria, C.A.R. Sa de Melo, Phys. 

Rev. B 55, 15153 (1997) 
* R. Shankar, Rev. Mod. Phys. 66, 129 (1994) 
^ G. Eilenberger, Z. Phys. 182, 427 (1965) 
° L.P. Gor'kov, T.K. Melik-Barkhudarov, Sov. Phys. JETP 

13, 1018 (1961) 



